﻿ Approximation of continuous-time systems with discrete-time systems The continuous-time systems are replaced by discrete-time systems even for the processing of continuous-time signals 1 Impulse invariance method 2 Step invariance method 3 Finite Difference Approximation (FDA) 4 Bilinear Transform 1 General method for quality evaluation of approximation for band-limited systems AD converter – Analog to Digital Converter DA converter –Digital to Analog Converter step, ω2T- sampling S– sampling frequency Hs- transfer function of the analog () a system to be approximated Hz- transfer function of equivalent () d digital system The digital system - minimizes error e Suppose the systemand input are band limited: Hjω= 0, for >ω () aM ω Xω= 0, for >ω () aM ω 3 System identification tδ -identifies a system, not () − necessarilly band limited sintπ () -identifies a band - limited tπ π with ω 0 ⇒exterior of unit disc |z|=r>1 14 segment [-π/T, π/T) on imaginary axis ⇒ one wrapping on the unit disc • In order to avoid aliasing errors in the frequency response of the digital system obtained, the frequency response of the analog system should be completely included in frequency bandwidth -π/T , π/T • Band-limited analog system T/π = 0 of 1/ 2,center 1/ 2, 0zσ> () disc σ= ⇒ = = Imaginary axis , 0 1/ 2,center 1/ 2, 0sj z=ω () 11jT+ω zxjy== + 2 = 3511jTT−ω +ω () • When an analog stable system is approximated, ideally the equivalent digital system should also be stable • Analog system and digital system have the same frequency response if the imaginary axis on s-plane becomes the unit circle on the z-plane •Not true for this method 36 π Ωπ Ω≤ or 36 22 ⎛⎞ 211⎛⎞ += and =1 very closexy z⇒− ⎜⎟ ⎜⎟ 22⎝⎠ ⎝⎠ ππ ≤⇔ω=Ω≅Ω≤ Hω limited at : tg () aM MMMTωΩ 36 36 ⋅ ω 36 2 very high value!!!!!⇒ω sM≥ ω ⇒ ω= ω Sampling theorem 2 at limit 36′′ sM s sω≥ Very high sampling frequency ⇒appropriate for modeling low frequency systems (electromechanical systems) 37 Example: first order LPF 1 = time constant: RCτ= 0ω 1 Hs= () a 1s+τ T T+τ () dHz= τ 1− 1z− T+τ 38 Finite difference approximation 39 Finite difference approximation 40 Finite difference approximation higher accuracy ! 41 4 Bilinear Transform t 1dy =ττ⇒ = () () ( ) ( ) axt yt x d H s=⇒ ∫ dt s −∞ Input signal Area A(ABCD)~ integral I: numerical methods (trapezoidal rule) nn nT ττ () nIxd= ∫ 1nT− 42() nT 1IynTynT xd=−−=ττ () ( )() () n∫ 1nT− () ⎤ +− 1xnT x n T TAB CD AD⎡ () ( ) () + () ⎣⎦ ≅= 22 T −= +− 11yn yn xn xn⇒− [] [ ][] [ ] () 2 1− 1Tz+ == ()() 1Hz Hs− 121daz− − = − 121sz− 1Tz+ • Bilinear transform: 1− 21z− == () () dasHzHs−1; 121z− − 1sTz=+ + 4311Tz− Relation of the s and z planes for the bilinear transformation method T + 11s− 21z− 2 = sz=⇒ 1− 1TTz+ 1s− 2 22 TT⎛⎞ ⎛⎞ 1+ σ+ ω ⎜⎟⎜⎟ 22⎝⎠ ⎝⎠ z=22 TT⎛⎞ ⎛⎞ σ+ ω 1− ⎜⎟⎜⎟ 22⎝⎠ ⎝⎠ ⇒ 01 zσ> ⇒ 0 1, right half plane (s plane) outside the unit disczσ> ⎩ 45 Relation between frequency responses for the bilinear transform Analog and digital frequencies connection: 2Ω tgω= or 2T Tω 2arctgΩ= 2 46 Tω 1Tjω+ 2arctgj 22jΩ 1sj z e re=ω⇒ = =⋅ = Tω 1j− 2 Tω rz== Ω=1 and 2 arctg 2 Imaginary axis (s⇒ plane) unit circle (z plane) Analog and digital frequencies: 2TΩω Ω= or: 2 arctg tgω= 22T 47 Distorted digital frequency response due to non-linear relation 48 between frequencies